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ABSTRACT 



Context. The mean UV extinction law for the Small Magellanic Cloud (SMC) is usually taken as a template for low-metallicity 
galaxies. However, its current derivation is based on only five stars, thus placing doubts on its universality. An increase in the number of 
targets with measured extinction laws in the SMC is necessary to determine its possible dependence on parameters such as metallicity 
and star-forming activity. 

Aims. To measure the UV extinction law for several stars in the quiescent molecular cloud SMC B 1-1. 

Methods. We obtained HST/STIS slitless UV spectroscopy of a 25"x25" field of view and we combined it with ground-based NIR 
and visible photometry of the stars in the field. The results were processed using the Bayesian photometric package CHORIZOS to 
derive the visible-NIR extinction values for each star. The unextinguished Spectral Energy Distributions (SEDs) obtained in this way 
were then used to derive the UV extinction law for the four most extinguished stars. We also recalculated the visible-NIR extinction 
for the five SMC stars with preexisting UV extinction laws. 

Results. The UV extinction law for four SMC Bl-1 stars within several pc of each other differs significantly from star to star. The 
2175 A bump is moderately strong in one, weak in two, and absent in the fourth. 

Key words. Dust, extinction — Galaxies: Magellanic Clouds — Ultraviolet: ISM — Stars: early-type 



1. Introduction 

In the 1980s IUE was extensively used to derive the UV extinc- 
tion law for a large number of stars in the Milky Way and the 
Magellanic Clouds. It became clear then that th ere were large 
variat ions among sightlines in the three galaxies. ICardelli et al.l 
(1989) used several tens of stars to derive a one-parameter fam- 
ily of extinction laws for the Galaxy. The number of useful tar- 
gets was even lower for the LMC and even mor e so for the 
SMC. The often quoted SMC extinction law of iPrevot et aT] 
(11984 was based in only 3 stars (AzV 398, AzV 18/Sk 13, and 
Sk 191). Sk 191 turned out to be a poor choice since it is actually 
an unredd ened object, leaving only t wo useful objects. The re- 
analysis of Gordon & Clavton ( 1998) a dded two further t argets, 
AzV 214 and AzV 456/Sk 143, and iGordon etail d200l com- 
pleted the current list of SMC stars with a derived UV extinction 
law with a fifth target, AzV 23/Sk 17. 

The characteristic that is most often quoted about the SMC 
extinction laws is t he absence of the 2175 A w ide absorption 
bump, which in the Fitzpat rick & Massal ([l990) parameteriza- 
tion translates into a zero or very small C3 term. However, this 
is not true of all SMC extinction laws since AzV 456/Sk 143 
clearly shows the absorption feature. The different location of 
AzV 456/Sk 143 in the SMC "wing" as c ompared with that 
of the other stars in the SMC bar prompted IGordon & Clavtod 
(1998) to compute two different SMC extinction laws, one for 
the wing and one for the bar, based on one an d thref] stars, re- 
spectively. The results of IGordon et al.l (120031) also suggest an- 
other possible difference between the two regions: Their mea- 



1 As we have already mentioned, Gordon et al. (2003) added another 
star, AzV 23/Sk 17, also located in the bar. 



sured Ry values for the four bar stars are in the typical range of 
2.40-3.30 while that of the wing star is, 2.05 + 0.17. The exis- 
tence of such detected variations in a small sample is an obvious 
sign of the need for a better sample to quantify the prevalence 
of the different extinction laws and to understand the underlying 
reason for such variations. Another reason why a better sample 
is needed is the relatively low valu es of E(B - V) for the five 
stars, 0-15-0.26 (IGordon et al.ll2003l) . which make the measured 
Ry values have large uncertainties. 

Why should we care about the SMC extinction law? Besides 
the importance that it has by itself, because the SMC is our best 
template for low-metallicity galaxies. Indeed, it was assumed for 
some time that the absence of the 2175 A bump in the SMC ex- 
tinction law was a metallicity effect but the detection of such 
a feature in Azv 456/Sk 143 seems to invalidate that hypothe- 
sis. Several explanations are possible, such as the existence of 
metallicity gradients in the ISM or the effect of "star-formation 
activity", a term that refers to the possible destruction of dust 
grains by UV radiation or shocks from massive stars. Until we 
have a better sample we will likely not be able to decide which 
explanation is the most plausible one. 

We selected as our target SMC Bl-1, a cold molecular cloud 
with no signs of star formation located at the southe rn end of the 
SMC bar. It has a molecular mass of 2.4- 10 4 M dRubio et al.l 
2004) as derived from the CO(1-0) emission line and a virial 
mass determination using observations done with the Swedish 
ESO Submilllimeter Telescope (SEST) at La Silla Observatory. 
Of the stars with existing UV extinction laws, the closest one to 
SMC Bl-1 is AzV 18/Sk 13, located -14' away. 
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2. Observations and data description 

2.1. STISdata 

SMC Bl-1 was observed with the Space Telescope Imaging 
Spectrograph (STIS) aboard the Hubble Space Telescope (HST) 
on 27 May 2004. We obtained four imaging exposures with 
STIS/NUV-MAMA, two imaging exposures with STIS/FUV- 
MAMA, and one slitless spectral exposure with STIS/NUV- 
MAMA. Both MAMA detectors have no read noise and, in 
imaging mode, cover a field of 25"x25" with 1024x1024 pix- 
els. The NUV imaging exposures used two different filters, 
F25CN270 (centered at 2702 A, see Fig.Q) and F25CN182 (cen- 
tered at 1895 A) with total exposures times of 100 s for each 
filter. The FUV imaging exposures used a single filter, F25QTZ 
(centered at 1559 A) with a total exposure time of 200 s. The 
slitless spectral exposures, which cover an area slightly bigger 
than that of the images, were obtained with the NUV objective 
prism (Fig. [2j with a total exposure time of 916 s. The NUV ob- 
jective prism has a variable dispersion that goes from «5 A/pixel 
at 1750 A to a45 A/pixel at 3000 A dSTSrill2007l) . Our spectral 
exposures were obtained in the shadowed part of the orbit and 
with the F25SRF2 filter in order to minimize the effect of the 
geocoronal background. 

There are 15 stars that are easily detected in all three UV 
imaging filters and in the spectral exposures. They have been 
labelled in Figs. Q] and [2] in order of ascending y coordinate. 
Note that the stars have not been placed on an absolute reference 
frame, so it is possible that the RA+dec grid is offset by as 
typical with coordinates derived from the header information in 
HST images. 

2.2. NIR imaging 

Deep NIR JHK S imaging was obtained with ISAAC at the 8.2 m 
VLT telescope at Paranal Observatory on 24-25 September 1999 
(63.C-0329(A)). We used the short wavelength arm equipped 
with a 1024 x 1024 1024 pixels Hawaii Rockwell array, with 
a spatial resolution of (X'148"/pix and a total FOV of 2;5x2.'5. 

The observations were done in a series of 6 frames (NDIT), 
each individual frame with a 10 sec integration time (DIT) in 
each filter. The individual frames were coadded and a 60 sec 
image was stored. A 10 position dither mosaic was done with 
a separation of 15" using the same observing strategy to cover 
the SMC Bl-1 region. In each filter, for every 10 minutes of on- 
source imaging we interleaved sky frame observations. These 
sky frames were chosen in a field with faint stars and no extended 
emission located 300" south of the position of SMC Bl-1. The 
sky field was observed in a similar way as the source frames. 
The procedure was repeated until we achieved a final integration 
time of 3600 s in J, H, and K s . 

To produce the final images, each image was dark corrected, 
flat fielded, and sky subtracted, and then median averaged and 
combined using IRAF procedures. The final images were regis- 
tered with respect to the J image by means of several common 
stars. The final ISAAC/VLT images cover a 2!5x5'0 area and 
have limiting magnitudes of about / = 21, H = 21, and K s = 22, 
respectively. 

We identified the 15 stars observed in the HST field in 
the VLT/NIR images and performed aperture photometry using 
the IRAF/DAOPHOT package. We used the standard stars ob- 
served on both nights to calibrate the photometry in K s (observed 
the first night) and in H and J (observed the second night). 
Additionally, we selected several stars from the 2MASS IR im- 



Table 1. ISAAC/VLT NIR photometry. 



Star 


J 


H 


K s 


01 


19.414 ± 0.069 


19.402 ± 0.072 


19.525 ± 0.094 


02 


20.002 ± 0.072 


20.090 ± 0.074 


19.968 ±0.104 


03 


20.086 ± 0.073 


20.006 ± 0.076 


19.963 ±0.100 


04 


19.864 ± 0.038 


19.797 ± 0.045 


19.647 ± 0.041 


05 


18.836 ± 0.066 


18.740 ± 0.068 


18.718 ±0.082 


06+07 


17.987 ± 0.065 


17.994 ± 0.068 


18.000 ± 0.079 


08 


18.436 ± 0.066 


18.317 ±0.068 


18.235 ± 0.080 


09 


18.424 ± 0.066 


18.378 ± 0.068 


18.307 ± 0.081 


10 


18.812 ±0.066 


18.759 ±0.068 


18.739 ± 0.082 


11 


18.464 ± 0.067 


18.400 ± 0.069 


18.426 ± 0.083 


12 


20.246 ± 0.055 


19.382 ± 0.035 


18.877 ± 0.041 


13 


19.402 ± 0.024 


19.400 ± 0.034 


19.258 ± 0.045 


14 


17.658 ± 0.065 


17.165 ±0.067 


17.066 ± 0.078 


15 


17.236 ± 0.065 


16.681 ±0.067 


16.566 ± 0.077 



ages with good photometry (AAA label) and compared their cat- 
alog magnitudes with our JHK S photometry. The agreement be- 
tween our photometry and 2MASS is good to 0.01 magnitudes. 
The NIR photometry in given in Table Q] Two of the stars (06 
and 07) are too close together in the ground-based data to be 
resolved, so unresolved magnitudes are provided. 

2.3. Visible photometry 

The Johnson-Cousins UBVI photometry of the sources presen t 
in the STIS expo sures was obtained fromlZaritskv et alj {2002) 
using the VizieR dOchsenbein et alj2000l) J/AJ/123/855 cat alog. 
13 of the 15 stars are resolved in the|Zaritskv et al. (2002) data 
while the pair 06+07 appears as a single unresolved source. 

3. Processing of the SMC B1-1 data 

Our goal in this paper is to obtain the UV extinction law of as 
many of the stars in the SMC Bl-1 region for which we can 
measure their spectra as possible. The strategy to accomplish 
that is: 

1. Measure the effective temperature (7W) and optical/NIR ex- 
tinction of the 14 sources which are point-like in the ground- 
based photometry. 

2. Extract the spectra for the 15 point sources in the NUV slit- 
less exposures. 

3. Select the sources with the highest optical/NIR extinction, 
since those objects are the ones for which the UV extinction 
should be more easily measured. 

4. Calculate the extinction law k(A) for those stars selected in 
the previous step. 

3. 1. Effective temperatures and optical color excesses 

The traditional mechanism to measure UV e xtinction laws is 
known as the pair method (Massa et al1 ll983l) and requires ob- 
taining spectra of two stars of the same spectral type, one with 
high extinction and the other one with low extinction. The fluxes 
of the two stars are corrected for distance effects if necessary, 
the low-extinction spectral energy distribution is divided by the 
high-extinction one, and the result is normalized by the extinc- 
tion difference between the two. The pair method requires the 
accurate measurement of the spectral type of the extinguished 
star, which may not be straightforward for dim objects, and the 
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Fig. 1. F25CN270 STIS NUV-MAMA image of the SMC Bl-1 region. The stars for which photometry was extracted are labelled. 



availability of data for a star with the same spectral type, which 
can lead to errors due to mismatches between stars. 

Alternatively, it is possible to measure extinction without re- 
ferring to a standard star by using a synthetic spectral energy 
distribution (SED) library and a numerical code that matches 
the available (spectro-)photometry to a grid with different stellar 
para meters (e.g. T e s, log and Z) and extinction am ounts and 
laws (Ma fz Apellanidl2004t iFitzpatrick & Massall2005l) . In this 
way, it is possible to derive the intrinsic or unextinguished SED 
directly from the data and, hence, calculate the extinction law 
directly by comparison with the observed spectral fluxes. Here 
we will follow t hat approach for the optical /NIR extinction us- 
ing CHORIZOS dMaiz Apellanizl2004l20"05ah . a Bayesian code 



that allows for the fitting of (spectro-)photometric data to an 
SED grid of up to five dimensions, including parameterized ex- 
tinction lawfl For the amount of extinction CHORIZOS uses as 
parameter £(4405 - 5495) = A(4405) - A(5495), the monochro- 
matic equivalent to E(B - V), because the former depends only 
on the amount and type of dust while the latter also depends on 
the spectral type of the sourcfl In other words, the same amount 
and type of dust in front of an O star and an M star will pro- 



2 CHORIZOS is a public sofware written in IDL that can be down- 
loaded from http://jmaiz.iaa.es 

3 Throughout this work, wavelengths will be expressed in A unless 
otherwise explicitly stated. 
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Fig. 2. STIS NUV-MAMA objective prism exposure of the SMC Bl-1 region. See Fig.Q]for orientation and field size. 



duce the same £(4405 - 5495) but different E(B - V), so using 
E(B - V) as an indicator of the amount of extinction can le ad to 
some counterintuitive results (see e.g. Ma ssev et al.l 120051) and 
should be avoided to characterize an extinction law. 

In order to use CHORIZOS to derive the intrinsic SED 
of the UV-bright objects in SMC Bl-1 it is useful to take 
into account different considerations. First of all, optical-NIR 
extinction laws have relatively simple functional forms and 
the Galactic laws can in principl e be accurately des cribed us- 
ing a single-family parameter (Car delli et all Il989lfl In this 



4 Recent works have determined that there may be more variations 
than previously expected in the NIR extinction law but the effect of 



wavelength range the published SMC extinction laws are not 
too differen t from the Galactic onefl: for example, the SMC- 
bar law of iGordon & Clavtonl ([1998) is quite similar to the 
ICardelli et all d 1989b law with #5495 = 2.8 for A > 3000 
A (i?54 95 = A(5495)/£'(44 05 - 5495) is the monoch romatic 
equivalent to Rv)- See also iMafz Apellaniz et al.l (120121) . In the 
UV, the extinction morpho logy is more complicated and, as 
iFitzpatrick & Massa (l2007h have shown, cannot be accurately 



those variations should be relatively small in our case since, as we will 
later see, our targets have all A K < 0.2 magnitudes. 

5 The outstanding exception being that of AzV 23/Sk 17, but see be- 
low. 
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Table 2. Second-run CHORIZOS results. Results for stars 14 
and 15 are not shown due to the poor quality of the fits. 



Star 


Teg 


£•(4405 - 5495) 


reduced x^- m 


01 


13 600 + 1300 K 


0.14 + 0.05 


1.1 


02 


12 600+ 900 K 


0.03 ± 0.03 


6.8 


03 


12 600+ HOOK 


0.06 ± 0.04 


1.8 


04 


13 500 + 3800 K 


0.19 + 0.15 


0.2 


05 


11900+ 800 K 


0.17 + 0.03 


3.7 


06+07 


13 700 + 700 K 


0.10 + 0.03 


3.8 


08 


23 500 ± 2900 K 


0.51 ±0.04 


1.5 


09 


16900 + 1300 K 


0.20 ± 0.03 


2.1 


10 


15 600 + 1300 K 


0.38 ± 0.03 


4.8 


11 


19 000+ 1900 K 


0.27 ± 0.04 


2.8 


12 


12400 + 3100 K 


0.14 + 0.04 


3.2 


13 


10400+ 800 K 


0.09 ± 0.06 


2.2 



described by a single-parameter family. Therefore, it makes 
sense to derive the intrinsic SEDs by excluding the photometric 
data with the shortest wavelengths because of the likely con- 
fusion between intrinsic (e.g. T^) and extinction-law effects. 
Thus, we will use the longest-wavelength UV photometric data 
(F25CN270) combined with the optical (UBVI) and NIR (JHK) 
data for that purpose. Second, given the location in the SMC 
and the observed magnitudes, we expect our stars to be early-to- 
mid B main-sequence stars with metallicity near log Z = -1.0. 
Also, given that the optical-NIR colors of non-supergiant B stars 
are relatively insensitive to gravity and metallicity, it will be ap- 
propriate to constrain our synthetic SED grid to main-sequence 
models with logZ = -1.0 and to use T eS , £(4405 - 5495), and 
extinction law type as the free parameters in the grid. Note that 
the best photometric discriminant for the effective temperature 
of B stars in the optical-NIR range is the Balmer jump and that 
with our filter selection we have two photometric points to its 
left (F25CN270 and U) and two points to its right (B and V), 
thus allowing us to measure temperature precisely. Finally, in 
order to use consistent zero-points for the different photometric 
systems involv ed, we will use the Joh nson sensitivity curves as 
determined by Mafz Apellaniz (2006) and the zero-point system 
of lMafz Apellaniz! (120071) . 

Taking into account the considerations in the previous para- 
graph, we executed an initial CHORIZOS run for each of the 14 
point-like visible sources using as input the F25CN270 + UBVI 
+ JHK photometry and as SEDs Kurucz atm ospheres with 
main-s equence gravities and log Z = - 1 .0 and ICardelli et al.l 
( 1989) extinction laws, leaving three free parameters to fit (r e ff, 
£(4405 - 5495), and R 54 g 5 ). The results of that initial run show 
that all 14 sources have £(4405 - 5495) lower than 0.60 and 
only three have values above 0.25. Given the low extinction of 
the 1 1 sources with £(4405 - 5495) < 0.25, their values of #5495 
are mostly unconstrained, as expected in such a case when the 
spectral types are not known a priori. For the three more extin- 
guished stars we derive values of ^5495 between 2.9 and 3.5 with 
uncertainties of ss 0.3 i.e. similar to those of the average Galactic 
and SMC extinction laws. This result prompted us to do a sec- 
ond and final CHORIZOS run with the same photometric data 
an d parameters but using a fixed extinction law, the SMC-bar 
of iGordon & Clavtonl (Q998). The obtained values for T e g and 
£(4405 - 5495) in the final CHORIZOS run are similar to those 
of the first one but with slightly lower error bars. The results of 
the second run are shown in Tableland plotted (for the four stars 
with the largest mean values of £(4405 - 5495)) in Figs.|3]and|4] 
As expected, all stars have temperatures consistent with being B 



starfl Their extinction-corrected luminosities are all consistent 
with being main-sequence stars. 

3.2. Spectral extraction 

In order to extract the multiple spectra prese nt in the objec- 
tive pr ism data (Fig.EJl, we used MULTISPEC (Mafz Apellaniz 
2005b), a software package specifically designed for the ex- 
traction of multiple spectra from slitless HST exposures of 
crowded fieldfl MULTISPEC works by fitting multiple spa- 
tial profiles for each column in the spectral exposure. We used 
v2.0 of the code, which allows for the u se of tabulated spa- 
tial profiles and which has been applied by Knig ge et al.l (2008) 
to extract UV spectra from STIS/FUV-MAMA G140L expo- 
sures of 47 Tuc. For the flux c alibration of the data we used 
iMafz Apellaniz & Bohlinl d2005l) and included time-dependent 
sensitivity corrections. 

The extracted spectra for stars 08, 09, 10, and 1 1 are shown 
i n Fig. |H along with the mod e CHORIZOS SED, which assumes 
a IGordon & Clavtonl d!998l) SMC-bar extinction law. The ex- 
tracted spectra have been binned in wavelength in order to obtain 
a uniform S/N in each bin. The observed and model spectra have 
relatively similar values of f v in their common wavelength range, 
indicating that the extinction law f or those stars canno t be to o 
different from the SMC-bar one of IGordon & Clavtonl £1998). 
Note, however, that there are some readily apparent differences: 
stars 08 and 11 show a weak 2175 A absorption structure and 
the measured UV fluxes for stars 09 and 10 are larger than the 
model ones. 



3.3. Source selection 

As previously mentioned only three point sources (08, 10, and 
12) have £(4405 - 5495) > 0.25 and a fourth one (09) has a 
slightly lower value. Those will be the four sources for which 
we will calculate the UV extinction law since for the rest the 
relative uncertainties on k(A) are too large to yield a useful result. 
Note that previous studies have derived extinction laws for stars 
with even lower values of £(4405 - 5495), but that was possible 
because of the availability of data with higher S/N that what we 
are using here. 

3.4. UV extinction-law calculation 

Finally, we obtain the UV extinction law k(A) which, following 
the usual practice, we normalize as: 



k(A) = fc(5495) + 



E(A - 5495) 
£(4405 - 5495) 



(1) 



and then consider only the second term, which is for 
A = 5495 A and 1 for A = 4405 A. Note, that, as previously 
explained, in this paper we use monochromatic rather than filter- 
integrated values for the definition of the extinction law. The 
mode value of k(A) - £(5495) is obtained by dividing the ex- 
tracted spectrum for each star by the unextinguished version of 
the CHORIZOS -derived mode SED and normalizing the result. 



6 Very recently, Mafz Apellaniz et al. (2012) have derived a new op- 
tical/NIR extinction law for 30 Doradus. Using that extinction law does 
not introduce fundamental changes to the results shown here. 

7 MULTISPEC is a public sofware written in IDL that can be down- 
loaded from http://jmaiz.iaa.es 
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Fig. 3. Contour plots showing the likelihood in the r e fj-£'(4405 - 5495) plane for the four stars in the second CHORIZOS run with 
the highest mean value for £(4405 - 5495). The contour intervals are linearly spaced with the minimum and maximum at 0.05 and 
0.95 times the peak (mode) value, respectively. The star shows the location of the mode (not the mean) and the plots on each axis 
show the integrated likelihood for that parameter. 



A similar procedure can be done for the filter-integrated quan- 
tities derived from the measured F25QTZ and F25CN182. Both 
spectral and photometric results are shown in Fig. [5] One readily 
apparent result is that there is a decrease in the intensity of the 
2175 A absorption structure as we follow the stellar sequence 1 1 
-» 08 -> 09 -> 10. 

One practical way of characterizing the UV extinction 
law is by using the six- parameter (ci, q, C3, C4, xq, y) 
Fitzpatrick & Massa dl990h recipe, which gives: 



k(A) - £(5495) = c\ + c 2 x + aDix, x Q , y) + caF(x), 



(2) 



where x = l/A (usually expressed in fi '), 



D(x, xq, 7) = 



(A' 2 



+ X 



2 y 2' 



and 



F(x) 



0.5392(x - 5.9) 2 + 0.05644(jc - 5.9) 3 , x > 5.9 
0, x < 5.9. 



(3) 



(4) 



We wrote a y 2 - r ninim ization code in order to fit such a 
Fitzpatrick & Massa] d 19901) six-parameter function to the spec- 
tral (1700-2900 A) and photometric (F25QTZ and F25CN182) 
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Fig. 4. MULTISPEC and fixed extinction-law CHORIZOS results for the four stars with the highest £(4405 - 5495). The red and 
blue data points with error bars correspond to the observed UV-optical-NIR photometric magnitudes, with the vertical error bars 
showing the uncertainties and the horizontal ones the approximate extent of the filters. Only the blue data poin ts points were used 
for the CHORIZOS fits. The continuous FUV-to-NIR line shows the CHORIZOS mode SED, which assumes a lGordon & Clayton! 
(1998) SMC -bar extinction law throughout the full wavelength range, while the NUV-only histogram shows the MULTISPEC 
extraction. The green stars show the synthetic photometry that corresponds to the CHORIZOS mode SED. The CHORIZOS values 
for r e jf and £(4405 - 5495) are given for each star. 



data for each of the four stars with higher extinction. When doing 
so, one has to be careful to account not only for random errors 
arising from the finite S/N of the data but also for systematic 
ones arising from possible spectral mismatches (in our case, this 
corresponds to the uncertainty in the r e jf and £(4405 - 5495) 
values derived by CHORIZOS). The results are shown and plot- 
ted in Fig. [5] 

As previously noted, there is a sequence of intensities of the 
2175 A absorption structure (measured by C3) in the order 1 1 — » 



08 — » 09 — > 10. It should be also pointed out that the values of 
the parameter that measure the FUV slope (04) are anomalously 
high. However, those values are not very significative (their er- 
rors are very large) because they are based mostly on a single 
photometric point (F25QTZ) centered at relatively long FUV 
wavelengths. 

Note that for our calculations we do not perform a prior sub- 
traction of the Galactic contribution to the extinction. As de- 
scribed by IZagurvl (120071) . for the Magellanic Clouds that can- 
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not be done from the data itself unless one assumes a smooth 
spatial distribution of the Galactic cirri, something that the MIR 
observations do not support. Therefore, the UV extinction laws 
derived here are in principle a combination of a Galactic and 
an SMC components. However, the existence of stars with very 
low extinction in the vicinity [e.g. star 02 has £(4405 - 5495) = 
0.03 + 0.03] argues in favor of the dominance of the SMC com- 
ponent. An additional point against the importance of a Galactic 
dust component is the wide range in reddenings in our sample, 
with values from 0.03 to 0.51 mag. Note that at a distance of 
e.g. 1 kpc to a hypothetical Galactic cloud, 10" correspond to 
just 10000 A.U. To our knowledge, no Galactic cirrus shows 
such large variations on such small scales. The low value of 
£(4405 - 5495) for star 02 is also in agr eement with the fore- 
ground reddening of 0.037 meas ured by ISchlegel et al] (l 998) 
and w ith the more recent work of lSubramaniam & Subramanianl 
(l20Toh . who find that most of the SMC has a Galactic foreground 
extinction lower than Ay — 0.1 mag. 

4. Looking back at the five SMC stars with 
preexisting extinction laws 

In order to check that our methods are consistent with pre- 
vious works, we used CHORIZOS to calculate the values of 
£(4405 - 5495) and R 5495 for the five SMC stars with preexist- 
ing ex tinction laws. For that purpose, we used the lGordon et alj 
(2003) UBVRI and the 2MASS JHK S photometry for those 
starqj. In this case we used a newer version of CHORIZOS that 
allows the stellar parameters to be T e $, luminosity class, and 
metallicity. By luminosity class, we mean a quantity that runs 
from 0.0 (hypergiants) to 5.5 (ZAMS) and, for a given T e ff and 
metallicity, yields a luminosity similar to the one derived from 
spectral classification (with Ia+ translate d as 0.0 and Vz as 5.5) . 
The S EDs in this case were TLUSTY dLanz & Hubenvll2003l 
|2007|) . Three parameters were left fixed: (from the spec- 
tral type), metallicity (SMC), and the logarithm of the distance 
(4.782). Three parameters were allowed to vary: Luminosity 
class, ^495, and £(4405 - 5495). For AzV 398 we excluded 
the K s photometry from the fit after an initial run because the 
star appears to have an excess in th at band compared t o the 
TLUSTY models, likely a wind effect. iCardeili et ail (fl989b ex- 
tinction laws were used. Results are shown in Table [3] Note that 
the values there for ^5495 and £(4405 - 5495) are monochro- 
matic quantities while Ay is filter-integrated. 

The reduced x^ a m all cases is quite good, indicating the 
good quality of the photometry. The (photometry-derived) lumi- 
nosity classes are all between 1.0 and 2.0, as expected for su- 
pergiants. Our results for £(4405 - 5495) and ^5495 differ sig- 
nificantly fro m the band-integrate d ones [E(B - V) and Ry, re- 
spectively] of iGordon et ail ([2003). The explanation for the dif- 
ferences in £(4405 - 5495) is likely to be caused by the fact that 
we are measuring absolute extinctions while they are measur- 
ing extinctions relative to a comparison star. Indeed, our values 
for £(4405 - 5495) are larger than their values of E(B - V) by 
0.01 to 0.10 magnitudes, as expected under such circumstances. 
However, we also measure values of ^5495 that are systemati- 
cally larger by ~0.6. That is harder to assign to the inclusion 
of an additional ^5495 = 3.1 component in our values (the one 
responsible for the additional extinction), since that would pull 
them towards 3.1 and not increase all of them independently of 



Note that, being supergiants, they are much brighter than the SMC 
Bl-1 stars. Hence, they all have 2MASS photometry with small uncer- 
tainties. 



^5495 . The difference may be attributed to the different method- 
ologies and to the fact that our quantities are monochromatic and 
theirs are band-integrated. 

The uncertainties in Table [3] (and in Table [2] as well) are 
formal (or random) ones. There are also systematic uncer- 
tainties caused by the mismatch between the SEDs used and 
the re al ones and due to the problems with the Cardell i et ail 
(119891) extinction laws (fundamentally, the use of a seventh 
deg ree polynomial as a functional form in the optical range , 
see [M afz Apellani zetail 120071: iMafz Apellaniz & Sotal 120081: 
IMafz Apellaniz et~ 20121) . However, an analysis of the pho- 
tometric residuals indicates that the Balmer jump seems to be 
well fit. Furthermore, these objec ts have £(4405 - 54 95) < 0.5, 
so the systematic effects of the Carde llTet alJ (Q989) extinction 
laws should be smalfl Therefore, we do not expect the system- 
atic uncertainties to be larger than the random ones, especially 
for R 5495 . 

The weighted mean of the ^5495 values for the four SMC bar 
stars is 3.36 + 0.07 and the results are consistent with the four 
real ^5495 values being the same (three cases within one sigma, 
one between one and two sigma). Therefore, the results for those 
four stars and the three more extinguished stars in SMC Bl-1 
are consistent with a visible-NIR extinction law not too different 
from the typical Galactic one, which has ^5495 = 3.1 - 3.2. On 
the other hand, the result for AzV 456/Sk 143 does indeed seem 
different, since its ^5495 is ~5 sigmas lower. 



5. Discussion 

In this paper we have calculated for the first time the UV ex- 
tinction law for four stars in SMC Bl-1, a quiescent cloud in 
the Small Magellanic Cloud. For those stars and for the other 
five stars with previously existing UV extinction laws we have 
also calculated their £(4405 - 5495) and ^5495 from visible-NIR 
photometry. 

For the four stars in SMC B 1 - 1 . all located within a few pc of 
each other, we have found significant variations in the extinction 
law and values. The latter is simply a manifestation of the small- 
scale structure in the dust distribution while the former implies a 
different dust composition. The most clear variations in the UV 
extinction law are in the strength of the 2175 A bump, which 
goes from non-existing (as it is the case for the four of the SMC 
bar stars with preexisting laws) to a case of moderate strength by 
Galactic standards though still weaker than the previously stud- 
ied SMC bar star (c 3 = 1.76 + 0.30 for the new target, 2.57 ±0.22 
for the old one). C3 does not appear positively correlated with 
£(4405 - 5495) (the most and the least extinguished stars have 
weak bumps while the two stars with intermediate extinctions 
represent the extremes in bump strength) which argues against a 
foreground Galactic cloud as the source of the extinction. 

Despite the significant UV variations, the SMC visible-NIR 
extinction law appears to be more uniform, at least within 
the current measurement errors. The variations in ^5495 have 
been attributed to differen t dust grain sizes for a long time 
dBaade & Minkowskilll937l) . All observed eight SMC bar stars 
with well-measured ^5495 have values compatible with the typ- 
ical Galactic one and only the SMC star measured in the wing 
has a somewhat lower R5495. However, there are no clear cases 
of high values (^5495 ~ 5), such as thos e found in some Galactic 
Hn regions (see e.g. lArias et alJ [2006). That is not surprising, 
given that the Galactic cases with ^5495 are found within the 



9 A quite different situation c ompared to the more e xtinguished case 
of Pismis 24-1 case discussed in Mafz Apellaniz et al. (2007). 
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c 2 = 
c 3 = 

Q = 

x = 
Y = 



Star 1 1 
1.80± 0.20 (r.)± 0.11 (s.) 
1.44± 0.06 (r.)± 0.09 (s.) 
1.76 ± 0.30 (r.)± 0.14 (s.) 
9.26 ± 2.36 (r.)± 0.95 (s.) 
4.76 ± 0.02 (r.) ± 0.01 (s.) 
0.85 ± 0.07 (r.)± 0.01 (s.) 




Fig. 5. Observed ext i nction laws (histogram used for the spectral data, points with error bars for the photometric data) and fitted 
Fitzp atrick & Massal (1 19901) extinction laws (continuous line) for each of the four stars in Fig. |4] The values of the fitted parameters 
are given in each case with their random and systematic errors indicated by "r." and "s.", respectively. 



bright areas of their H n r egi ons and that is not the case for any 
of the analyzed SMC star{3 

The existence of the 2175 A bump in stars located in rela- 
tively quiescent regions of the SMC argues in favor of its absence 
being caused by the presence of intense UV radia tion fields and 
stellar winds, which wo uld destroy its carrier (see lClavton et al.l 
2003: ISofia et all 12006). Nevertheless, our results also suggest 
that the final explanation may be more complex. Otherwise, the 
four stars in SMC Bl-1 would show the bump in their extinc- 
tion laws. Also, since extinction appears to be patchy and as- 
sociated with local dust clouds, cj, should be correlated with 



AzV 214 is in NGC 346 but outside the central bright nebulosity. 



£(4405 - 5495). Neither of those two effects appear in on our 
results. 
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Table 3. CHORIZOS results for the five stars with preexisting UV extinction laws. 



Star 






luminosity class 


R 


5495 


£(440f 




- 5495) 


A, 




red. y L - 


AzV 18/Sk 13 


19000 


K 


1.19 


±0.01 


3.93 


±0.64 


0.178 


± 


0.022 


0.696 ± 


0.040 


1.9 


AzV 23/Sk 17 


18 000 


K 


1.01 


±0.01 


3.41 


±0.08 


0.253 


± 


0.003 


0.878 ± 


0.015 


1.9 


AzV 214 


23 000 


K 


1.80 


±0.03 


3.09 


±0.20 


0.246 


± 


0.007 


0.772 ± 


0.034 


2.2 


AzV 398 


28 000 


K 


1.93 


±0.03 


3.57 


±0.28 


0.322 


± 


0.018 


1.162 ± 


0.036 


1.3 


AzV 456/Sk 143 


28 000 


K 


1.19 


±0.02 


2.63 


±0.14 


0.337 


± 


0.011 


0.901 ± 


0.027 


0.6 
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